Formation of a Partially- Screened Inner Acceleration Region in 
Radio Pulsars: Drifting Subpulses and Thermal X-Ray Emission 

from Polar Cap Surface 

Janusz Gil 1 ' 2 , George Melikidze 1,3 & Bing Zhang 2 
ABSTRACT 

The subpulse drifting phenomenon in pulsar radio emission is considered 
within the partially screened inner gap model, in which the sub-Goldreich- 
Julian thermionic flow of iron ions or electrons coexists with the spark-associated 
electron-positron plasma flow. We derive a simple formula that relates the ther- 
mal X-ray luminosity L x from the spark-heated polar cap and the E x B subpulse 
periodicity P 3 (polar cap carousel time). For PSRs B0943+10 and B1133+16, 
the only two pulsars for which both P 3 and L x are known observationally, this 
formula holds well. For a few other pulsars, for which only one quantity is mea- 
sured observationally, we predict the value of the other quantity and propose 
relevant observations that can confirm or discard the model. Then we further 
study the detailed physical conditions that allow such partially screened inner 
gap to form. By means of the condition T c /T s > 1 (where T c is the critical tem- 
perature above which the surface delivers a thermal flow to adequately supply 
the corotation charge density, and T s is the actual surface temperature), it is 
found that a partially-screened gap (PSG) can be formed given that the near 
surface magnetic fields are very strong and curved. We consider both curvature 
radiation (CR) and resonant inverse Compton scattering (ICS) to produce seed 
photons for pair production, and find that the former is the main agency to 
produce gamma-rays to discharge PSG. 
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1. Introduction 

Pulsar radio emission typically occurs in a form of periodic series of narrow bursts of 
radiation. These burst often have a complex structure, in which higher order periodicities 
can be found. The individual pulses consist of one, few to several subpulses. In some pulsars 
the subpulses demonstrate a very systematic drift across the pulse window. If the pulses are 
folded with the basic pulsar period then the drifting subpulses form amazingly spectacular 
patterns called drift-bands. 

The phenomenon of drifting subpulses is a long standing puzzle in the pulsar research 
and its solution would likely result in deeper understanding to the nature of pulsar radiation. 
It is generally believed that this phenomenon is inherently associated with the so-called 
inner acceleration region above the polar cap, in which the magnetospheres plasma does not 
corotate with the neutron star surface. The first model based on this idea was proposed 
by Ruderman & Sutherland (1975, RS75 henceforth). The predictions of RS75 model were 
successfully compared with a handful of pulsars known to show this phenomenon at that 
time. A decade later, Rankin (1986, R86 hereafter) compiled a list of about 40 drifting 
pulsars, but drifting subpulses were still regarded as some kind of exceptional phenomenon. 
However, recently Weltevrede, Edwards, & Stappers (2006, WES06 hereafter) presented the 
results of a systematic, unbiased search for subpulse modulation in 187 pulsars and found 
that the fraction of pulsars showing drifting subpulse phenomenon is likely to be larger 
than 55%. They identified 102 pulsars with drifting subpulses in their sample, with a large 
fraction of newly discovered drifters. The authors concluded that the conditions required 
for the drifting mechanism to work cannot be very different from the emission mechanism 
of radio pulsars. WES06 then suggest that the subpulse drifting phenomenon is an intrinsic 
property of the pulsar emission mechanism, although drifting could in some cases be very 
difficult or even impossible to detect due to insufficient signal-to-noise ratio. It is therefore 
essential to attempt to unravel the physical conditions that can lead to formation of an inner 
acceleration region above the polar cap that could lead to development of the subpulse drift 
phenomenon. 

The classical vacuum gap model of RS75, in which spark-associated sub-beams of sub- 
pulse emission circulate around the magnetic axis due to E x B drift of spark plasma fila- 
ments, provides the most natural and plausible explanation of drifting subpulse phenomenon. 
However, despite its popularity, it suffers from the so-called binding energy problem. Namely, 
under canonical conditions the surface charges (ions or electrons) are likely to be directly 
pulled out of the surface so that a pure vacuum gap is difficult to form. The alterna- 
tive steady flow polar cap models, the so-called space-charge-limited flow models (Arons & 
Scharlemann 1979; Harding & Muslimov 1998), cannot give rise to the intermittent "spark- 
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ing" behavior, which seems necessary to explain the subpulse drift phenomenon in radio 
pulsars. Gil & Mitra (2001, GM01 hereafter) revisited the binding energy problem of RS75 
model and argued that the formation of the vacuum gap (VG) is, in principle possible, al- 
though it requires a very strong non-dipolar surface magnetic fields, much stronger than a 
canonical dipolar component inferred from the observed spindown rate. Once the binding 
is strong enough to prevent the thermionic emission at the full space charge limited flow, 
the inevitable E x B drift of plasma filaments will result in the observable subpulse drift 
phenomenon. It has been known for a long time that in order to allow all radio pulsars 
to produce electron-positron pairs (the necessary condition for coherent radio emission), 
the near-surface magnetic fields must include multipole components dominating over global 
dipole field (Ruderman & Sutherland 1975; Arons & Scharlemann 1979; Zhang, Harding, & 
Muslimov 2000). Growing observational evidence of non-dipolar structure of surface mag- 
netic field 1 accumulates, and the suggestion that such a sunspot-like fields form during the 
early proto-neutron star stage has been proposed 2 (e.g. Urpin & Gil 2004). Gil & Melikidze 
(2002, GM02 hereafter) calculated the non-dipolar VG model for 42 drifting subpulse pul- 
sars tabulated by R86 and argued that VG can be formed in all considered pulsars, provided 
that the actual surface magnetic field was close to 10 13 G independently of the value of the 
canonical dipolar magnetic field. 

Although the binding energy problem could be, at least in principle, resolved by as- 
suming an appropriately strong surface magnetic field, yet another difficulty of the RS75 
model was that it predicted a much too fast E x B drift rate. Motivated by this issue, 
Gil, Melikidze, & Geppert (2003, GMG03 hereafter) developed further the idea of the inner 
acceleration region above the polar cap by including the partial screening by a sub-Goldreich- 
Julian thermal flow from the surface due to the spark-associated polar cap heating. This 
idea was first introduced by (Cheng & Ruderman 1980, CR80 henceforth), who argued that 
even with thermionic ions included in the flow, the condition above the polar cap is close 
to that of a pure vacuum gap. A similar model was also invoked by Usov & Melrose (1995, 
1996). GMG03 reanalyzed this model and argued that a thermostatic self-regulation should 
keep the surface temperature just few percent below the critical ion temperature at which 
the gap potential drop is completely screened. This results in more than 90 % of screening 



1 It is worth mentioning that although RS75 implicitly assumed non-dipolar surface magnetic fields to 
treat the j-B pair production processes, in their calculations of many other physical quantities (such as the 
surface charge density) they still used dipolar form, presumably for the sake of simplicity. 

2 Anothcr possible mechanism of creating small scale anomalies of surface magnetic fields was proposed 
by Geppert, Reinhardt & Gil (2003). They argued that due to a Hall drift instability, the poloidal magnetic 
field structures can be generated from strong subsurface toroidal fields. 
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due to thermionic emission. Given a similar gap height, the actual potential drop, and hence 
the E x B drift rate, is about 10 % of that of the pure vacuum case. This is still above the 
threshold for the magnetic pair production. Similar to the avalanche pair production cascade 
introduced by RS75, the discharge of the growing potential drop above the polar cap would 
also occur in the form of a number of sparks. In this paper we call such an inner accelerator 
a partially screened gap (PSG henceforth). The latest XMM- Newton observation of the 
drifting pulsar PSR B0943+10 (Zhang, Sanwal, & Pavlov 2005, ZSP05 hereafter) reveals a 
possible hot spot with the surface area much smaller than the conventional polar cap. This 
is consistent with the polar cap heating from such a PSG, which at the same time gives the 
right E x B drift rate. This lends strong support to the PSG model. 

In this paper we study PSG model in greater detail and explore the physical conditions 
for the model to work. We also apply this model to a new set of 102 pulsars from WES05 
and show that it can work in every case, provided that the surface non-dipolar magnetic 
field is strong enough, even stronger than 10 13 G suggested by GM02 and GMG03. Our new 
treatment is a combination of those used in GM02 and GMG03. Our working hypothesis 
is that drifting subpulses manifest the existence of a thin inner acceleration region, with an 
acceleration length scale much shorter than the polar cap size. The ultra-high accelerating 
potential drop discharges via a number of localized E x B drifting sparks, as has been 
envisioned by RS75 . These sparks produce isolated columns of electron-positron plasma that 
stream into the magnetosphere to generate radio-beams of the observed subpulse emission 
due to some kind of plasma instability (see the §8 for more discussion). Due to charge 
depletion with respect to the co- rotational Goldreich- Julian (1969) value, sparks experience 
an unavoidable E x B drift with respect to the polar cap surface. As a consequence, the 
spark-associated sub-beams of radio emission perform a slow circumferential motion that 
is responsible for the observed subpulse drift. This model, which is often called "a pulsar 
carousel model" (Deshpande & Rankin 1999), is examined in this paper 3 . We are particularly 
interested in the thermal effect associated with the surface bombardment by back-flowing 
particles produced by sparks. The intrinsic drift rate (manifested by the tertriary drift 
periodicity) and the polar cap heating rate (manifested by the thermal X-ray luminosity) 
should be correlated with each other, since they are determined by the same value of the 
accelerating potential drop. The properties of drifting subpulses are discussed in §2 and the 
properties of charge depleted acceleration region above the polar cap are discussed §3. We 



3 Other suggestions of subpulse drifting as phenomena occurring outside the inner gap have been made 
(Kazbegi, Machabeli & Melikidze 1991; Spitkovsky & Arons 2002; Wright 2003; Gogoberidze, Machabeli, 
Melrose, et al. 2005; Fung, Khechinashvili & Kuijpers 2006), but the connection between the radio drifting 
rate and the X-ray properties in those models is not yet clear, and we do not discuss them in the current 
paper. 
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find a specific relationship between the appropriate observables and conclud that it holds 
for a number of pulsars for which good quality data is available (especially PSR B0943+10, 
ZSP05). It turns out that this relationship depends only on the observational quantities and 
thus it is a powerful tool for testing this theoretical model. Although the number of pulsars 
that have all necessary data for such testing is small at the moment, the clean prediction 
holds the promise to ultimately confirm (or discard) the PSG model in the future. In §§4-7 
we analyze this model in a more detailed manner to investigate the microscopic conditions 
(e.g. near-surface configuration, radiation mechanism, etc) that are needed to form such 
PSGs. The paper is summarized in §8. 



2. Drifting subpulse phenomenon 

Let us briefly review the characteristic properties of the drifting subpulse phenomenon, 
in which subpulses typically change phase from one pulse to another in a very organized 
manner, to the extent that they form apparent driftbands of duration from several to a few 
tens of consecutive pulses. Usually more than one drift band appears, and the separation 
between them measured in pulsar periods ranges from about 1 to about 15 (Backer 1976; 
Rankin 1986). The subpulse intensity is systematically modulated along drift bands, decreas- 
ing towards the edges of the pulse window. In some pulsars, however, only periodic intensity 
modulations are observed, without any systematic phase change. These pulsars were iden- 
tified as those in which the line-of-sight cuts through a beam centrally (Backer 1976), thus 
showing a steep gradient of the polarization angle curve (e.g. Lyne & Manchester 1988). On 
the other hand, the clear subpulse driftbands are typically found in pulsars associated with 
the line-of-sight grazing the beam, thus showing a relatively flat position angle curve (Backer 
1976; Rankin 1986; Lyne & Manchester 1988). It is then obvious that both these kinds of 
subpulse behavior, i.e. systematic drift in phase and phase stationary intensity modulation, 
represent the same phenomenon, namely, rotation of beams of subpulse emission around the 
pulsar magnetic axis. The picture of rotating subbeams is also supported by spectral prop- 
erties of the observed radiation. The separation between drifting subpulses in a given pulse 
show dependence on frequency similar to those of profile components and/or overall pulse 
width, while the observed periodicities related to patterns of drifting subpulses are indepen- 
dent of radio frequency. This seems to exclude all frequency-dependent plasma effects as a 
plausible source of the drifting subpulse phenomena. Also, a similarity of drift patterns at 
different radio frequencies strongly suggest that the radiation of drifting subpulses is related 
to a relatively stable system of isolated plasma filaments originating near the polar cap. 

The observational characteristics of drifting subpulses described briefly above suggest 
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then strongly an interpretation of this phenomenon as a number of isolated sub-beams of 
radio emission, spaced more or less uniformly in the magnetic azimuth, and rotating slowly 
around the magnetic axis. In terms of RS75 model this means that the subpulse-associated 
sparks perform a circumferential motion around the pulsar beam axis. The most spectacular 
confirmation of this "carousel model" was presented by Deshpande & Rankin (1999, 2001) 
and Asgekar & Deshpande (2005), who performed a sophisticated fluctuation spectra analysis 
of single pulse data from PSR B0943+10 and detected clear spectral features corresponding 
to the rotational behavior of subpulse beams. These results indicate clearly that sparks do 
not vanish after passing through the pulse window in the form of drifting subpulses, but 
they continue to drift around the polar cap and reappear in the same phase of the pulse 
window after a definite and measurable time period. The rotational features appearing in 
the fluctuation spectra of PSR B0943+10 have extremely high values of Q ~ 500 (defined 
as the central frequency divided by the width) confirming a circulational nature of the 
phenomenon. Similar feature with Q ~ 100 was found by Asgekar & Deshpande (2005) in 
PSR B0834+06. Quite recently Gupta, Gil, Kijak, et al. (2004) revealed a system of subpulse 
beams circulating around the polar cap in PSR B0826-34. The high values of Q suggests 
strongly that the trajectories of the spark motion associated with drifting subpulses should 
be closed. It is natural to conclude that the circulating sparks would form on average in a 
ring on the polar cap corresponding to a conal beam of radio emission originating at some 
height above the polar cap (Gil, Kijak & Seiradakis 1993; Kijak, & Gil 1998). 

The "pulsar carousel model" based on E x B drifting sparks is widely accepted by radio- 
observers but theorists usually treat it with a great deal of reservation. Besides the binding 
energy problem (see §4), they usually point out problems with spark stability and mem- 
ory, which is indeed necessary to explain the systematic subpulse drift within the carousel 
model. This is a very important but difficult problem that requires sophisticated computer 
simulations. Here we would only raise some phenomenological arguments. Following RS75, 
Gil & Sendyk (2000; GS00 henceforth) argued that both the characteristic spark dimension 
and the typical distance between the adjacent sparks are approximately equal to the gap 
height. At any instant the polar cap should be populated by as many sparks as possible. 
Indeed, spark discharges should develop wherever the accelerating potential drop exceeds the 
threshold for the magnetic pair production. Such a maximum packing should stabilize an 
instantaneous arrangement of sparks on the polar cap surface. Each spark develops and dies 
over a time scale of microseconds, so the question is what agency makes them to reappear in 
the same place (such reappearance is necessary to explain the regular subpulse drifting). It 
seems that a natural reason is thermal emission of charges at the base of each spark is heated 
by the return bombardment. Moreover, some electrons and positrons from the outflowing 
spark plasma column would still be near the polar cap surface when the returning potential 
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drop again exceeds the pair formation threshold. These charged particles would initiate the 
very next discharge (see Asseo & Melikidze 1998; Melikidze, Gil & Pataraya 2000, and GS00 
for more detailed discussion). Therefore, the only net motion of sparks is the slow E x B 
drift across the planes of the surface magnetic field. For the line-of-sight grazing the pulsar 
beam (corresponding to sparks operating at the polar cap boundary), this circumferential 
drift motion will be observed as a systematic subpulse drift, while for a more central line- 
of-sight trajectories only longitude-stationary modulation of the subpulse amplitude can be 
observed. 



3. Charge depleted inner acceleration region 

The inner acceleration region above the polar cap may be a result of deviation of the 
local charge density p from the co-rotational charge density (Goldreich & Julian 1969, GJ69 
hereafter) 

f2 ■ B s B s 

PGJ = ~2^- * ^P' (1) 
where the positive/negative sign corresponds to ^Fe ions/electrons. These two possibilities 
correspond to antiparallel and parallel relative orientation of the magnetic and spin axes, 
respectively. As mentioned before, there exists growing evidence (see Urpin & Gil 2004, 
and references therein) that the actual surface magnetic field B s is highly non-dipolar. Its 
magnitude can be described in the form 

B s = bB d (2) 

(Gil & Sendyk 2000, GS00 hereafter), where the enhancement coefficient b > 1 and 

B d = 2 x 10 12 (PP_ 15 ) 1/2 G (3) 

is the canonical, star centered dipolar magnetic field, P is the pulsar period in seconds and 
P_i 5 = P/10~ 15 is the period derivative. For the purpose of further considerations in this 
paper we introduce here a convenient measure of the surface magnetic field 

|^ = 0.046(PP_ 15 ) a5 & (4) 

Pq 

where B q = 4.414 x 10 13 G is the so called quantum magnetic field (Erber 1966). 

The polar cap is defined as the locus of the so-called open magnetic field lines that 
penetrate the light cylinder (GJ69). Generally, the polar cap radius can be written as 

r p = 1.45 x io 4 p-°- 5 r a5 cm (5) 
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(GSOO), where the factor 6~ 5 describes squeezing of the polar cap area due to the magnetic 
flux conservation in non-dipolar surface fields (as compared with the GJ69 dipolar configu- 
ration). One should realize, however, that the polar cap radius r p expressed by the above 
equation is only a characteristic dimension. In fact, in the presence of strong non-dipolar 
surface magnetic field the actual polar cap can be quite irregular in shape. Also, the light- 
cylinder radius and thus the canonical polar cap radius depends on the unknown particle 
inertia (e.g. Michel 1973). Therefore, r p as expressed in equation (5) can be used as an order 
of magnitude estimate of the actual polar cap radius. 

The observationally deduced polar cap radii (see Section 3.2) are often much smaller than 
the canonical GJ69 value. This seems consistent with strong non-dipolar surface magnetic 
field, that is large b in equation (5). Given a strong enough near-surface magnetic field (the 
required condition will be discussed below in §§4-7), charge depletion in the acceleration 
region above the polar cap can result from binding of the positive l®Fe ions (at least partially) 
in the neutron star surface. Positive charges then cannot be supplied at the rate that 
would compensate the inertial outflow through the light cylinder. As a result, a significant 
part of the unipolar potential drop (GJ69, RS75) develops above the polar cap, which can 
accelerate charged particles to relativistic energies and power the pulsar radiation. The 
characteristic height h of such an acceleration region is determined by the mean free path of 
pair-producing high energy photons. In other words, the growth of the accelerating potential 
drop is limited by the cascading production of an electron-positron plasma (e.g. RS75, 
CR80). The accelerated positrons would leave the acceleration region, while the electrons 
would bombard the polar cap surface, causing a thermal ejection of ions, which are otherwise 
more likely bound in the surface in the absence of additional heating. This thermal ejection 
would cause partial screening of the acceleration potential drop AV corresponding to a 
shielding factor 

rj = 1 - A (6) 
Pgj 

where p\ is thermonically ejected charge density (see also eq.[18]) and 

= v^B s h 2 , (7) 

is the accelerating potential drop 4 , where B s is defined by equation (2) and h is the model 
dependent height of the acceleration region (see section 5). The above expression for AV 
is the solution of Poisson equation for a thin PSG, assuming that dr]/dh = within the 
acceleration region. 



4 In this paper we ignore slight reduction of AV due to general relativistic effect of frame dragging, which 
can affect the potential drop by as much as 27 percent (see Zhang et al. 2000; Gil & Melikidze 2002). 
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The cascading sparking discharge in PSG works analogously to the scenario envisioned 
by RS75. The only difference is that the accelerating potential drop is much smaller due to 
the partial screening by thermionic ions with charge density p x . Additional positive charge 
density p + is supplied by the spark produced positrons. The potential drop is completely 
screened when the total charge density p — p\ + p + reaches the corotational value (eq.[l]). 
Then the spark plasma escapes from the region as discussed by Asseo & Melikidze (1998). 
However, during a short escaping timescale h/c ~ 1CT 8 s the thermionic ions still maintain 
the partial screening as determined by equation (6). Thus, the returning potential drop 
quickly grows to the level defined by equation (7) and the spark begins to develop again. 
This results in an intrinsically intermittent nature of PSG and an unsteady particle flow 
into the magnetosphere, which has an important physical consequences for generation of the 
coherent pulsar radio emission (see section 8). 

In the stationary observer's frame the co-rotational charge density pgj = — • B/2-7TC 
is associated with the co- rotational electric field E cor (r) = —(17 x r) x B(r)/c and the 
magnetosphere co-rotates with the velocity v cor = c(E cor x T5)/B 2 , where B is the magnetic 
field (dipolar) and ft is the pulsar spin axis (Q = 2n/P). Now let us introduce a thin, charge 
depleted region extending to about h above the polar cap surface, with the accelerating 
potential drop described by equation (7). Since the charge density within this region is lower 
than the co-rotating value pcj, an additional electric field AE appears in this region, which 
makes the spark plasma to drift with a velocity v d = c/B 2 (AE x B), as observed in the 
co-rotating frame. To estimate the values of AE let us follow the original method of RS75 
and consider the closed contour abfea within the gap, as marked in Figure 4 of RS75. Since 
the potential drop along segments ab and bf vanishes (AV a b = AVbf = 0), one then has 
AV ae = AVf e , where AV ae corresponds to the tangent electric field AE, while AV e / is just 
the acceleration potential drop expressed by equation (7). However, we can consider another 
segment f'e' (parallel to fe) and use the same argument to demonstrate that AV ae > = AV ae 
for any arbitrary pair of points e and e'. This means that AE ~ within the RS75 gap, 
except at the boundary region within about h from the polar cap boundary, where h is the 
height of the gap. The perpendicular electric field AE grows rapidly from essentially zero 
to the corotational field E cor in the boundary region. The situation can be slightly different 
if the gap height h (length of segments fe and f'e') varies between the pole and the polar 
cap edge, although in the thin gap case of RS75 (h <C r p ) such variations cannot be large 
(Ah/h <C 1). In this case AVf e ^ AVf> e > and AV ae ^ AV ae i and some residual AE can exist 
over the entire polar cap area. In any case, the tangent electric field is strong only at the polar 
cap boundary where AE = 0.5AV/h = rjiji / cP)B s h (see Appendix A in GMG03 for details). 
Here B s is the surface magnetic field (eqs. [3] and [2]), AV^ is the accelerating potential drop 
(eq. [7]), h is the gap height (eqs. [19] and [20]), and rj is the shielding factor (eqs. [6], [9] and 
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[16]). This electric field causes that discharge plasma at the polar cap boundary performs 
a slow circumferential motion with velocity i><j = cAE/B s = i]7rh/P. The time interval to 
make one full revolution around the polar cap boundary is P 3 ~ 27ir p /v d . One then obtains 
the so called tertiary drift periodicity 



^3 r p 



P 277/1' 



(8) 



which is the proper measure of the intrinsic drift rate, at least at the boundary of the polar 
cap. If the plasma above the polar cap is fragmented into filaments (sparks) which determine 
the intensity structure of the instantaneous pulsar radio beam, then, at least in principle, 
the tertiary periodicity P 3 can be measured/estimated from the pattern of the observed 
drifting subpulses (e.g. Deshpande & Rankin 1999; Gil & Sendyk 2003). In practice P 3 is 
very difficult to measure (mainly because of aliasing which is a severe problem even in pulsars 
with high signal-to-noise ratio, e.g. DR99, GS03) and its value is known at the moment only 
in few cases. It is much easier to measure the primary drift periodicity P 3 , which in high 
signal-to-noise ratio is just a distance between the observed drift bands measured in pulsar 
periods P (there are also clever techniques that allow to measure P 3 even in cases with very 
low signal-to-noise ratio; see Edwards & Stappers 2002, WES06). Since P 3 = NP 3 (e.g. 
RS75), where N « 2ixr p /2h is the number of sparks contributing to the drifting subpulse 
phenomenon observed in a given pulsar (GS00), then one can write the shielding factor in 
the form 

1 P 

which depends only on relatively easy-to-measure primary drift periodicity P 3 (WES06). 
This equation can be compared with the definition of shielding factor (eq. [6]) in order to 
derive the amount of thermally ejected iron ions or electrons (depending on the sign of charge 
of the polar cap). Apparently, the shielding parameter rj should be much smaller than unity. 
Note that in the RS75 pure VG model one has 77 = 1, which implies the predicted subpulse 
drift at least an order of magnitude too fast as compared with observations. This means 
that the conditions within the natural inner acceleration region above the polar cap should 
greatly differ from the pure vacuum gap proposed by RS75. 



3.1. Thermostatic self-regulation of the potential drop 

GMG03 argued that because of the exponential sensitivity of the accelerating potential 
drop AV to the surface temperature T s , the actual potential drop should be thermostatically 
regulated (see also CR80). In fact, when Ay is large enough to ignite the cascading pair 



- 11 - 



production, the back-flowing relativistic charges will deposit their kinetic energy in the polar 
cap surface and heat it at a predictable rate. This heating will induce thermionic emission 
from the surface, which will in turn decrease the potential drop that caused the thermionic 
emission in the first place. As a result of these two oppositely directed tendencies, the quasi- 
equilibrium state should be established, in which heating due to electron bombardment is 
balanced by cooling due to thermal radiation. This should occur at a temperature T s slightly 
lower than the critical temperature above which the polar cap surface delivers thermionic 
flow at the corotational (GJ69) charge density level (see GMG03 for more details). 

The quasi-equilibrium condition is <tT s 4 = / ym c c 3 n, where 7 = eAV/m e c 2 is the Lorentz 
factor and AV is the accelerating potential drop (eq. [7]). Here 

n = n G j - rii = r]n GJ (10) 

is the charge number density of back-flowing particles that actually heat the polar cap surface, 
i] is the shielding factor (eq. [6]), n; is the charge number density of thermally ejected flow 
and 

n G j = Pcj/e = 1.4 x lO 11 ^ ^- - 5 cm" 3 (11) 

(eq. [1]) is the corotational (GJ69) charge number density 5 . It is straightforward to obtain 
an expression for the quasi-equilibrium polar cap surface temperature in the form 

T s = (2 x 10 6 K) p-°- 25 P _^ v - 5 b - 5 hl- 5 , (12) 

where /13 = h/10 3 cm is the normalized height of the acceleration region. This height is 
model dependent and we discuss two possible models below (eqs. [21] or [22]). Growing 
evidence suggests that in actual pulsars T s is few MK (see Table 1). Thus, from the above 
equation one can infer that b ~ 77 -1 , where the shielding parameter should be much lower 
than unity (see discussion below eq. [7]). As a consequence, the enhancement coefficient 
b = B s /B d (eq. [2]) should be much larger than unity (at least of the order of 10). This 
is very consistent with the binding energy problem discussed in Section 4. Indeed, the 
extremely high surface magnetic field helps to resolve this problem. 

3.2. Interrelationship between radio and X-ray signatures of drifting subpulses 

The predicted thermal X-ray emission luminosity from the polar cap with temperature 
T s is L x = aT^A bo \, where A ho \ = nr 2 and r p is the actual polar cap radius (eq. [5]) and 



5 Note that in the original RS75 paper the charge number density corresponds to purely dipolar magnetic 
field (b = 1), while in our approach the surface magnetic field is highly non-dipolar (b ^> 1). 
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a is the Stefan-Boltzmann constant. Thus L x = 1.2 x 10 32 (P-15 / P 3 )(i]h/r p ) 2 erg/s, which 
can be compared with the spin-down power E = IQti = 3.95/45 x 10 31 P_is/P 3 er g/ s > where 
/ = /4 5 10 45 g cm 2 is the neutron star moment of inertia (in what follows we assume that 
J 4 5 = 1. We can now use the equation (8) in the form f]h/r p m O.5P/P3 and derive the 
thermal X-ray luminosity from the polar cap heated by sparks, i.e. 

^ = 2.5xl0^^j , (13) 

where P3 is the tertiary periodicity in drifting subpulses pattern (pulsar carousel time), which 
is equal to the time interval needed for the discharge plasma to make one full revolution 
around the perimeter of the polar cap. 

One can also derive the X-ray luminosity efficiency (with respect to the spin down 
luminosity) 




(14) 

We can see that equations (13) and (14) depend only on the observational data of the radio 
pulsars. It is particularly interesting and important that both equations above do not depend 
on details of the sparking gap model (i],b,h). Thus, we have found a simple relationship 
between the properties of drifting subpulses observed in radio band and the characteristics of 
X-ray thermal emission from the polar cap heated by sparks associated with these subpulses. 
For PSR B0943+10, which is the only pulsar for which both P 3 = 37P, (DR99) and L x = 
5 x 10 28 erg s _1 ~ 5 x 10~ 4 P, (ZSP05) are measured, the above equations hold very well. 
In few other cases for which the circulational periodicity is measured 14 < P 3 /P < 37 (see 
Table 1), which gives L x /E ~ (llo.s) x 10~ 3 - Such correlation between the X-ray and the 
spin-down luminosities is a well-known and intriguing property of rotation powered pulsars 
(Becker & Triimper 1997; Possenti, Cerratto, Colpi, et al. 2002). It has been suggested that 
this correlation is a characteristic of a magnetospheric radiation (e.g. Cheng, Gil, & Zhang 
1998; Zhang & Harding 2000). Here we suggest that it can also be a characteristic property 
of the polar cap thermal radiation. Most likely, both mechanisms contribute at a comparable 
level to the observed X-ray luminosity. However, one should realize that if the estimates for 
PSRs 0809+74 and 0633+08 (Table 1) were correct, they would violate significantly the 
above correlation, so it should be treated with caution until confirmed or discarded. 



Using equations (8) and (12) we can write the polar cap temperature in the form 

-1/2 

P 



T s = (5.1 x 10 6 K)&V4pV4p-i/2 ( £| ) ^ (15) 
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where the enhancement coefficient b = B s / B d A pc /A ho \, A pc = nr^ and A ho \ = irr^ (see 
Section 2). Since Aboi can be determined from the black-body fit to the spectrum of the 
observed thermal X-ray emission, the above equations can also be regarded as independent of 
details of the sparking gap model and depending only on the combined radio and X-ray data, 
similarly as equations (13) and (14). This seems quite understandable, since it reflects the 
simple fact that within the sparking gap model both the intrinsic drift rate and the polar cap 
heating rate are determined by the same electric field. Therefore, the observational values 
of L x and P3/P should be related to each other through a combination of P and P that is 
proportional to the spin-down power E, as it is really so in equations (13) and (14). Let us 
summarize a set of assumptions that lead to these equations: (i) the electric field causing 
the drift is estimated at the polar cap boundary 6 (eq. [8]), (ii) the gap is partially shielded 
(Eq. [10]), and (iii) the surface magnetic field at the polar cap is highly non-dipolar (eq. [11]). 
If some or all of these assumptions were not valid, then the unknown parameters such as h, b 
and/or rj would appear in equations (13) and (14). Therefore, and observational verification 
of these equations seems to be an easy way to confirm or discard the PSG model of the inner 
accelerator in pulsars. 



3.3. Actual pulsars 

Table 1 presents the observational data and predicted values(computed from equations 
(13 - 15)) of a number of quantities for five pulsars, which we believe show clear evidence of 
thermal X-ray emission from spark heated polar caps. Besides the three cases with known 
P 3 (B0943+10, B0826-34 and B0834+06), we also included PSR Bl 133+16 (twin of PSR 
B0943+10) for which we estimated P 3 and argued that its value was actually measured but 
misinterpreted as the primary drift periodicity P 3 (shown in parenthesis in Table 1). PSR 
B0809+74, for which an estimate of P3 exits, is included. We also added the Geminga pulsar 
(B0633+17) for which thermal radiation from small polar cap was clearly detected. 

PSR B0943+10. This is the best studied drifting subpulse radio pulsars with P = 1.09 s, 
P_i 5 = 3.52, E = 10 32 erg s' 1 . DR99 clearly demonstrated that the observed subpulse drift 
is in this pulsar aliased and found the alias-resolved values of drift periodicities P3 = 37.4P 
and P3 = 1.86P, which gave the number of observed sparks N = P3/P3 = 20. In an 
attempt to detect thermal X-ray signatures of these sparks ZSP05 observed this pulsar 
with the XMM-Newton. They obtained a spectrum consistent with thermal BB fit with a 



6 This is different from the original RS75 model, where the electric field is estimated at the middle of the 
polar cap. Moreover, they used pure vacuum gap model. 
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bolometric luminosity L x rs 5 x 10 28 erg s" 1 and a polar cap surface area A ho \ = 10 7 (T s /3 x 
10 6 K)~ 4 cm 2 ~ (II04) x ^ 7 cm2 5 which was much smaller than the conventional polar cap 
area A pc = 6 x 10 8 cm 2 . This corresponds to the best fit temperature T s = 3.1 x 10 6 K (see 
Fig. 1 in ZSP05). The predicted values of L x and L x /E calculated from equation (13) and 
(14), respectively, agree very well with the observational data. The surface temperature T s 
calculated from equation (15) with b = A pc /Ab \ is also in very good agreement with the 
best fit. The shielding parameter rj « 0.1 is derived from equation (9). 

PSR B1133+16. This pulsar with P = 1.19 s, P_i 5 = 3.7, and E = 9 x 10 31 erg s" 1 is 
almost a twin of PSR B0943+10. KPG05 observed this pulsar with Chandra and obtained 
a spectrum consistent with thermal BB fit L x = 6.7 x 10 28 erg s -1 (we used their L x /E = 
3.6lg;? x 10-4 < cos^ with < cos9 >= 0.47), A hol = (0.5±g;|) x 10 7 cm 2 and T s « 2.8 x 
10 6 K. These values are also very close to those of PSR B0943+10, as should be expected for 
twins. Using equations (13) or (14) we can predict P3/P = 27^1 for B1133+16. Interestingly, 
Rankin (1986, R86 henceforth) gives the primary drift periodicity P3 = (5.3 ± 1.2)P for 
the conal components (probably aliased), but in the saddle between these components R86 
reports P3 = (25 ± 3)P. However, recently WES06, reported in this pulsar P3/P = 3 ± 2, as 
well as long period drift feature with (33 ± 3)P in the trailing profile component. We thus 
claim that this long period (shown in parenthesis in Table 1) is the actual tertiary periodicity 
P 3 /P = 33I3, which agrees very well with our prediction P3/P = 27I2. This interpretation 
is strongly supported by Fig. 10 in Nowakowski (1996), where the fluctuation spectrum of 
PSR B1133+16 shows prominent long period feature at 0.031 cycles/period, corresponding 
to 32 pulsar periods. It is worth noting that P3/P3 = 33 ± 3 is very close to 37.4 measured 
in the radio twin PSR B0943+10. Note also that the number of sparks predicted from our 
hypothesis is iV = P3/P3 = (33 ± 3)/(3 ± 2) = lllg 5 , thus it is quite possible that actually 
N is close to 20, as in the case of twin pulsar B0943+10. 

PSR B0826-34. With P = 1.85 s, P_ i5 = 1.0, E = 6 x 10 30 erg s _1 , this pulsar has 
P 3 « IP (highly aliased) and P 3 ss 14P, i.e. iV = 14 (Gupta, Gil, Kijak, et al. 2004; Esamdin, 
Lyne, Graham-Smith, et al. 2006). With these values we can predict from equation (13) that 
L x = 2 x 10 28 erg s^ 1 . This pulsar should be as bright as PSR B0943+10 and thus it is 
worth observing with XMM-Newton. The shielding parameter rj m 0.16 (eq. [9]). 

PSR B0834+06. With P = 1.27 s, P_i 5 = 6.8, and E = 1.3 x 10 32 erg s" 1 , this pulsar 
has P3 = 15P, and P3 = 2.16 P (Asgekar & Deshpande 2005), implying the number of 
sparks iV R3 15/2.16 ~ 7. From equation (13) we obtain L x = 37 x 10 28 erg s" 1 . This pulsar 
should be almost 8 times more luminous than PSR B0943+10, and we strongly recommend 
to observe it with XMM-Newton. The shielding parameter rj R3 0.07 (eq. [9]). 
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p^i? B0809+74- This is one of the most famous pulsars with drifting subpulses, with 
P = 1.29 s, -P-15 = 0.17, E — 3.1 x 10 30 erg s -1 and well determined value of the primary drift 
period P3 = 11 P. Recently, van Leeven, Stappers, Ramachandran et al. (2003) estimated 
the circulational tertiary period P3 > 150P. This is very long as compared with other pulsars 
considered above, but see the Geminga case below. Taking the lower limit P3 = 150 P we 
can predict from eq. (13) thermal X-ray luminosity at the level L x = 3 x 10 25 erg s _1 , which 
is a very low compared with other nearby pulsars. The efficiency ratio L x / E = 3 x 10~ 6 . 
The shielding factor rj = 0.014 (eq. [9]). 

PSR B0633+17 (Geminga). This pulsar with P = 0.237 s and P_i 5 = 10.97 is a 
problematic radio emitter, but well known for pulsating optical, X-ray and T-ray emission. 
Recently, Caraveo, De Luca, Mereghetti et al. (2004) detected thermal emission from an 
60 meter-radius polar cap heated to T s = (1.9 ± 0.3) x 10 6 K, with the bolometric thermal 
luminosity L x = 1.5 x 10 29 erg s _1 . Since the canonical radius of the polar cap is 300 m, 
then the amplification factor 6(300/60) 2 = 25 (eqs. [2] and [5]). Using the observed value of 
L x we obtain from equation (13) the predicted value of P3/P = 400, which is even longer 
than the minimum estimate in PSR B0809+74. Now taking this value and the estimate 
for b we obtain from equation (15) T s = 2.13 x 10 6 K, which is in very good agreement 
with observational range of the surface temperature (1.9 ± 0.3) x 10 6 K. The efficiency ratio 
L x /E = 4 x 10" 6 , close to the estimate in PSR B0809+74. 

Becker, Weisskopf, Tennant, et al. (2004) revealed the X-ray emission from a number 
of old rotation-powered pulsars. Among them is PSR B0823+26, which is well known to 
exhibit regularly drifting subpulses (see Table 1) with P 3 = (2.2±0.2)P (WES06). Such an 
old drifting pulsar is an ideal case to study the spark associated thermal radiation from hot 
polar cap. Becker, Weisskopf, Tennant, et al. (2004) found an upper limit of the thermal 
contribution from a hot polar cap as being T pc = 1.17 x 10 6 K, which seems low as compared 
with estimates given in Table 1. However, Becker, Weisskopf, Tennant, et al. (2004) used a 
full canonical polar cap with radius 299 meters to estimate T pc . In the actual non-dipolar 
surface magnetic field one should re-scale T s = 6 a25 T pc , where b = (r p /r pc ) 2 = B s /B^ 
(eqs. [2,5]). Since B& = 10 12 G in this case, b would be at least 10. This raises the lower 
limit of T pc to be above 2 x 10 6 K, consistent with our model. The same applies to another 
old pulsar B0950+08 studied by Becker, Weisskopf, Tennant, et al. (2004). 

4. Binding energy problem 

The phenomenon of drifting subpulses seems to indicate strongly a presence of the charge 
depleted acceleration region just above the polar cap. However, the existence of such region 
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strongly depends on the binding energy of |gFe ions, which are the main building blocks of 
neutron star surface (Usov & Melrose 1995; Lai 2001). If this cohesive energy is large enough 
to prevent thermionic emission, a charge depleted acceleration region can be formed above 
the polar cap. Normally, at the solid-vacuum interface, the charge density of outflowing ions 
is roughly comparable with density of the solid at the surface temperature kT s = e c , where 
e c is the cohesive (binding) energy and k = 8.6 x 10~ 8 keVK -1 is the Boltzman constant. 
However, in the case of pulsars, only the corotational charge density pgj can be reached, 
and the ^Fe ion number density corresponding to pgj is about exp(— 30) times lower than 
in the neutron star crust. Since the density of outflowing ions p\ decreases in proportion to 
exp(— e c /kT s ), one can then write Pi/pcj ~ exp(30 — e/kT s ). At the critical temperature 

T ' - m < 16 > 

the ion outflow reaches the maximum value p; = p GJ (Eq. [1]), and the numerical coefficient 
30 is determined from the tail of the exponential function with an accuracy of about 10%. 

Calculations of binding energies are difficult and uncertain (see Usov & Melrose 1995, 
1996 and Lai 2001 for detailed discussion). In this paper we use the results of (Jones 
1986, J86 henceforth), which were recommended by Lai (2001) in his review paper as 
more robust than others. J86 obtained e c =0.29, 0.60 and 0.92 keV for B s = 2, 5 and 10 x 
10 12 G, respectively. These values can be approximately represented by the function e c ~ 
(0.18keV)(-B s /10 12 ) ar G. Using equations (2), (3) and (16), this could be converted into the 
critical temperature 

7] ~ (0.7 x 10 5 K) f ^ (1-2 x 10 5 K) 6 7 (PP_ 15 ) » (lO 6 K) ( , 

(17) 

where we use the ratio B s /B q (eq. [4]) for the convenience of presenting the results in the 
graphical form (Fig. 1 and 2). Above this temperature the thermionic ion flow reaches the 
maximum GJ density at the surface, and the polar cap flow will be space charge limited. 
An acceleration potential can be still developed but the growth rate is slow compared with 
the PSG model we are developing in this paper. At temperatures below 7] charge-depletion 
would happen right above the surface, and an efficient acceleration region would form, which 
should be discharged in a quasi-steady manner by a number of sparks, similar to the original 
suggestion of RS75 for the pure vacuum gap case (see more detailed discussion in GMG03). 
The electron-positron plasma produced by sparking discharges co-exists with the thermally 
ejected ions, whose charge density can be characterized by the shielding factor (defined by 
eq. [6]) in the form 

7i' ~ 



77 = 1 — exp 
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As one can see from equations (6) and (18), at the temperature 7] = T s the shielding factor 
rj = (corresponding to fully developed space-charge limited flow with p x = pcj), but even a 
very small drop of T s below 7], much smaller than 10%, corresponds to creation of the pure 
vacuum gap with rj = 1 (pi = 0). Thus, the condition for partially screened charge depleted 
acceleration region can be written in the form T s < Tj, meaning that the actual surface 
temperature T s should be slightly lower (few percent) than the critical ion temperature 7], 
which for a given pulsar is determined purely by the surface magnetic field B s (eq. [17]). 
Practically, one can use 



to denote the condition of forming a PSG, with the caveat that in reality T s cannot be exactly 
equal to 7] but should be a few percent lower. 

As mentioned above, the polar cap surface can be, in principle, negatively charged 
(f2 • B > 0). In such a case (called "antipulsars" by RS75) the polar cap surface can deliver 
an electron flow. Following GMG03 we assume that this electron flow is also determined 
mainly by thermoemission, with the corresponding shielding factor r\ = 1 — p c / pqj = 1 _ 
exp[25(l — T c /T s )], where p c is the charge density of thermionic electrons. The critical 
electron temperature is 



(see GMG03 for more details), and in analogy to equation (19) the condition for creation of 
charge depleted acceleration region is T s = T e . Since the enhancement coefficient b » 1 (see 
the discussion following eq. [12]) so that B s / B q ~ 1 (eq. [4]), both 7] and T e can be close to 
T s ~ a few MK. Since 7] and T c are similar, we will only include the ion case (eq. [17]) in the 
following discussion, keeping in mind that our considerations are general and independent 
of the sign of the polar cap charge, at least qualitatively. 



So far, we have considered the macroscopic properties of the inner accelerator region in 
a way independent of details of the sparking gap model. Now we begin to discuss the micro- 
scopic properties of this model. In particular, we will present results of model calculations 
similar to those presented by GM02 but for a much larger sample of drifting subpulse pulsars 
and with the inclusion of partial screening due to thermionic emission from the polar cap 
surface. The condition for the formation of the gap used in GM02 was T c /T s > 1, while in 
this paper we used T c /T s = 1 (eq. [19]) as a result of the thermostatic regulation considered 
by GMG03. 



T s = T i 



(19) 




(20) 



5. PSG models of the inner acceleration region in pulsars 
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Growth of the accelerating potential drop (eq. [7]) is limited by the cascading production 
of relativistic electron-positron pair plasma in strong and curved magnetic field. RS75 derived 
a famous formula for their VG height assuming that the quasi-steady breakdown is driven by 
magnetic pair production induced by the curvature radiation seed photons. They used Erber 
(1966) approximation, which is valid in relatively low magnetic fields B s /B q < 0.1, which is 
not relevant for the magnetic field at the polar cap considered in this paper. In the strong 
surface magnetic field, i.e. B s > 5 ■ 10 12 G, the high energy photons with energy Ej = huo 
produce electron-positron pairs at or near the kinematic threshold hio = 2mc 2 / sin9, where 
sin^ = l P h is the photon mean free path for pair formation and 1Z = TZq ■ 10 6 cm 

is the radius of curvature of magnetic field lines. This regime is called the near threshold 
(NT) conditions (e.g. Daugherty & Harding 1983). Two VG models can be considered under 
the NT conditions: Curvature Radiation dominated Near Threshold Partially Screened Gap 
(CR-NTPSG) model and Inverse Compton Scattering dominated Near Threshold Partially 
Screened Gap (ICS-NTPSG) model, in which the potential drop is limited by pair production 
of the CR and the resonant ICS seed photons, respectively. Proper expressions corresponding 
to pure vacuum gap case have been derived by GM01 and GM02. Below we give generalized 
formulae by including partial screening due to the thermal flow from the polar cap surface 
(GMG03 and references therein). 

CR-NTPSG model 

In this model the cascading e _ e + pair plasma production is driven (or at least domi- 
nated) by the curvature radiation photons with typical energy hu) = (?>/2)h^c/TZ. where 
7 is the typical Lorentz factor of electrons/positrons moving relativistically along the local 
surface magnetic field lines with a radius of curvature 1Z. In the quasi-steady conditions the 
height h of the acceleration region is determined by the mean free path that is h ~ Z p h for 
pair production by energetic CR photon in the strong and curved magnetic field. Following 
GM01 (see also GM02) and including the partial screening effect (GMG03) we obtain 

h « h CR = (3.1 x 10 3 )nl /7 r)-V 7 b- 3 / 7 P 3 / 14 P:?l U cm. (21) 

ICS-NTPSG model 

In this model the cascading pair plasma production is driven (or at least dominated) 
by the resonant inverse Compton scattering (e.g. Zhang et al. 1997), with typical photon 
energy hu = 2 r yheB s /mc. In the quasi-steady conditions the height h of the acceleration 
region is determined by the condition h ~ l e where l c = 0.00277 2 ( J B S /10 12 G)- 1 (T S /10 6 K)- 1 
is the mean free path of the electron to emit this high energy ICS photon (ZHM00), -B s is 
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the surface magnetic field and T s is the actual surface temperature. Following GM01 and 
including partial screening (GMG03) we obtain 

h » h lcs = (5 x I^nfr^b^p-^PZU 2 cm. (22) 

In the above equations 7Z 6 = 71/ 10 6 cm is the normalized curvature radius of the surface 
magnetic field lines. 

Besides resonant ICS, also the non-resonant ICS with a characteristic energy ftw ~ 
max(7 2 /cT, •ymc 2 ) has been widely discussed in literature. This process usually dominates 
the gap dynamics in pulsars with moderate and low magnetic fields (Zhang et al. 1997; 
Hibschman & Arons 2001; Harding & Muslimov 2002). Since the polar caps in our model 
typically anchor a very strong curved magnetic field, the non-resonant ICS is unimportant 
and we ignore it in our detailed treatment. Throughout this paper ICS refers to resonant 
ICS, and the equations (21) and (22) will be referred as to CR and ICS cases, respectively. 

5.1. CR-NTPSG 

The results of model calculations for 102 pulsars with drifting subpulses (see Table 1) 
are presented in Figure 1 for the ion T\ critical temperature (eq. [17]). We have plotted B s /B q 
(left-hand side vertical axes) versus the pulsar number (which corresponds to a particular 
pulsar according to Table 1). The actual value of Bj B q (eq. [4]) for a given pulsar was 
computed from the condition T s = TJ, where T s is the actual surface temperature (eq. [12]), 
with the height h of CR driven acceleration region determined by equation (21). This condi- 
tion leads to the expression B s /B q = 126.4P- 2 °/ 29 ftg 0/2 y°/ 29 , which allows to sort pulsars 
according to decreasing value of pulsar period P, which is marked on the top of the Figure 1. 
The vertical axes on the right hand side are expressed in terms of the surface tempera- 
ture T s computed from equation (17). Different panels correspond to different normalized 
radii of curvature ranging from TZq = 0.1 to TZq = 0.009. Different symbols used to plot 
exemplary curves correspond to arbitrarily chosen values of shielding factor ranging from 
?] = 0.015 to rj = 0.15 (each curve represents the same shielding conditions rj for all 102 
pulsars considered). 

The vertical line in Figure 1 corresponds to PSR B0943+10 (Nl=41 in Table 1), which 
was observed using XMM-Newton by ZSP05. Three horizontal lines correspond to T s equal to 
about 2, 3 and 4 MK (from bottom to the top), respectively, calculated from equation (17). 
(More exactly, we marked T s = 7} = 2.08, 3.11 and 4.14 for B s /B q = 2.7, 4.8 and 7.2, 
respectively). Thus, the hatched areas encompassing these lines correspond to the range 
of surface temperatures T s = (3 ± 1) x 10 6 K deduced observationally for the polar cap of 
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B0943+10 from XMM-Newton observations (ZSP05, see their Fig. 1). 

5.2. ICS-NTPSG model 

Let us now consider the resonant inverse Compton scattering radiation model (ICS- 
NTPSG). The major difference between CR and ICS cases is the additional regulation of 
T s caused by the condition h ~ l e , where l e is inversely proportional to T s . In fact, the 
increase of the surface temperature causes the decrease of h, and hence, the decrease of 
AV. This makes the back-flow heating of the polar cap surface less intense. Let us es- 
timate the surface temperature in the ICS dominated case. The average Lorentz factor 
of electrons or positrons can be estimated by the gap height and the partially screened 
potential drop as 7 = (1.1 x 10 2 ) P l ^PZll 6 b 1 ^ 3 (T s /10 6 ) 2 / 3 i]^ 1 ^. Using the kinematic near 
threshold condition and the expression of the resonant ICS photons energy (see section 2) 
we get hw = (7.5 x 10~ 8 ) 7&(PP_i 5 ) -1 / 2 , which leads to another estimate for the aver- 
age Lorentz factor 7 = 2.5 (T s /10 6 )^ 3 TZ 1 / 3 - Combining these two expressions for 7 we 
find that for the case of the ICS dominated NTPSG the following relationship should hold 
7] = 8.14 x 10- 5 pV2p-V2 fo (t b /10 6 ) n 6 \ As a result, contrary to the CR case, the shielding 
factor i] is not a free parameter in the ICS case, and the analysis similar to that pre- 
sented in Figure 1 is not relevant. However, we can use the quasi-equilibrium condition 

<tT 4 = 'jm e c 3 n r], which leads to T s = (8.8 x 10 4 ) {^fqbP^^PZll 2 ^ . Using the expres- 
sions of i] and 7 we can then obtain 

T s = (1.5 x 10 4 K)7^" 1/4 p- 3 / 8 Pri 5 /8 6 3/4 . (23) 

Consequently, we can obtain the expression for the shielding factor 

77 = (1.2 x 10- 6 ) P^PZltn^'W^. (24) 

Thus, for a given pulsar (P, P) the values of T s and 77 are determined by the parameters of 
the local surface magnetic field TZq and b. In Figure 2 we plot T s as a function of B s / B q 
(eq. [4]) for different values of 1Z 6 and rj. The smooth lines in both panels correspond to 
the ion critical temperature, while other lines marked by different symbols shown in legends 
correspond to the actual surface temperature calculated for the CR case (upper panel) and 
the ICS case (lower panel). We look for partially screened solutions corresponding to the 
intersection of smooth and symbol-broken lines. 

As seen from the lower panel of Figure 2, equations (23) and (24) cannot be self- 
consistently satisfied. In fact, unlike in the CR case presented in the upper panel, in the 
ICS case the surface temperature T s ~ 10 6 K (eq. 23) is considerably lower than 7} (eq. 17). 
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Bringing T s close to 7] would imply 77 = 1 (pure vacuum), which is inconsistent with equation 
24. Thus, if the pair creation is controlled by ICS, then the acceleration region should be 
a pure vacuum gap rather than a PSG. This is consistent with previous results of GM02, 
who concluded that in the pure vacuum gap case ICS seed photons are more efficient in 
driving the cascading pair production. However, GM02 failed to understand that CR pair 
production mechanism produces more polar cap heating and thus makes the gap partially 
screened. In such a case the ICS pair production still occurs but it is no longer the dominant 
mechanism to control the properties (height, potential drop, etc) of the acceleration region 
above the polar cap. 



6. Thermal signatures of hot polar cap 

Knowing that ICS cannot contribute significantly to the pair production process in a 
PSG, let us consider polar cap heating by back-flow plasma particles created from CR seed 
photons. Using the near threshold condition in the form h = (4/3)(mclZ 2 /fry 3 ) = h 3 10 3 cm 
(see below eq. [4]) we can rewrite equation (12) in the form 

/D \ 2 / 7 

T s = (7.7 x 10 6 K)p- 2 /y/ 7 7^ /7 l^j . (25) 
Consequently, we can obtain efficiency of thermal X-ray emission in the form 

h = 0.15P 19 / 14 p-f ?? 8/7 7^ /7 (I 1 ) 1A , (26) 

where L x = aT s A^ Q \ is the X-ray luminosity from the hot polar cap with the surface area 
Awl = (eq. [5]) and E = IQQ = 3.95 x 10 31 P_i5/P 3 erg/s is the spin-down power. For 
the twin radio pulsars PSR B0943+10 and PSR B1133+16 the X-ray signatures of the hot 
polar cap are T s - 3 x 10 6 K, E « 10 32 erg/s and LJE ~5x 1(T 4 (see Table 1). Using these 
measurements we can infer from equations (25) and (26) an approximate condition in the 
form rfTZa ~ 6.5 x 1(T 5 . Since in PSR B0943+10 P 3 /P = 1.86 then according to equation 
(9) we obtain the normalized radius of curvature 1Z% = 0.009 in this pulsar (see next section). 
The microscopic equations (26) and (26) along with two macroscopic equations (14) and (15) 
provide full knowledge about the spark driven inner acceleration regions in pulsars. The PSG 
model behind these equations should be easily verifiable because they are either completely 
independent or weakly dependent on physical parameters of the acceleration region. 

Figures 3 and 4 present calculations of the efficiency L x /E and the temperature T s , 
respectively, for 102 pulsars from Table 1. One should emphasize that these calculations 
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utilized the two above equations and did not assume a priori an existence of the PSG in 
form of the condition T s = T c (eqs. [17] and [19]). Thus, these figures should be considered 
in association with Figure 1, which does take into account this condition. The vertical 
solid lines correspond to PSR B0943+10 for which T s — (3 ± 1) x 10 6 K (corresponding to 
B B /S q =4.8 in Fig. 1) and L x /E = (5 ± 2) x 1(T 4 (ZSP05). Due to weak dependence of the 
efficiency on the magnetic field we decided to use only one value of B s / B q =4.8 in Figure 3 to 
avoid overlapping of too many lines. The observational range of L x / E and T s are marked by 
the hatched belts in Figsure 3 and 4, respectively. Since L x /E in equation (26) depends on 
both P and P we have sorted pulsars in Figure 3 according to increasing value of E, which 
are marked on the top of the figure (N2 in Table 1). If L x /E ~ 10~ 3 1CT 4 for all pulsars, 
then one can say that the shielding factor rj should be larger for pulsars with larger E. Also, 
it follows from Figure 4 that if T s — (2 4) x 10 6 K in all pulsars, then the shielding factor 
7] should be larger for longer period pulsars. This is consistent with equation (9). 



7. Special case of PSR B0943+10 

The microscopic parameters of PSR B0943+10 are described in section 3.2 (with ap- 
propriate entry in Table 1). Let us now discuss this special case in more detail. PSR 
B0943+10 is the only pulsar for which both P 3 and L x are known observationally. As al- 
ready mentioned above, both the subpulse drift rate and the polar cap heating rate (due to 
subpulse-producing sparks) are related to the same value of the accelerating potential drop. 
In section 3, we found a "clear cut" formula involving two observables: the X-ray luminosity 
L x and the tertiary drift periodicity P 3 , which is independent of microscopic details of the 
acceleration region. This relationship holds very well for PSR B0943+10 (see Table 1). It is 
important to check whether this "model independent" approach is consistent with detailed 
model calculations, involving microscopic conditions prevailing in PSG. According to equa- 
tions (8) and (9) we have the shielding parameter rj = 0.09 and the complexity parameter 
r p/h = N/7T = 6.36, in consistency with GS00 (see also Gil, Melikidze, & Mitra 2002). Let 
us analyze in this respect our Figure 1, which allows to read off the physical parameters of 
the acceleration region for a given pulsar within a particular model rj. For PSR B0943+10, 
only the intersections of the vertical lines corresponding to Nl=41 (see Table 2) with the 
hatched belts, or even with the horizontal line T s = 3MK, are relevant. The above value of 
rj = 0.09 corresponds to the radius of curvature 1Z 6 ~ 6.5 x 10 5 ?7 2 ~ 0.009 (see the condi- 
tion below eq. [26]), and thus from the lower-right panel of Figure 1 we have B s / B q ~ 5 or 
B s ~ 2 x 10 14 G and T s ~ 3 x 10 6 K. This set of parameters can be also inferred from the 
upper panel in Figure 2. These values therefore are consistent with the appropriate entry for 
PSR B0943+10 in Table 1. This consistency between macro- and micro-scopic calculations 
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holds also for other pulsars from Table 1, as illustrated in Figure 3. This set of parameters 
can be also inferred from the upper panel in Figure 2. 

The surface temperature about 3MK (ZSP05), which implies B s = bB d ~ 2 x 10 14 G 
if one applies J86 theory (eq. [17]) to the case of B0943+10. The value exceeding 10 14 G 
may seem extremely high, but at least three radio pulsars have dipolar magnetic fields above 
10 14 G (McLaughlin et al. 2003, and references therein). The non-dipolar surface magnetic 
fields could be even stronger. The reason that the inferred non-dipolar field is even stronger 
than in the pure vacuum case found in GM02 for the ICS seed photons is that we are using 
a lower binding energy calculations of Jones (1986), which are more robust than those of 
Abrahams & Schapiro (1991) (see Lai 2001, for critical review). However, generally the 
surface magnetic field required to form CR-NTVG obtained by GM02 is still higher than 
that for CR-NTPSG obtained in this paper. 

In section 3.1 we assumed a pure black-body conditions at the polar cap. GM02 con- 
sidered deviations from the black-body conditions by introducing the so-called heat-flow 
coefficient k = Q ra d/(Qrad + Qcond), which described the amount of heat conducted beneath 
the polar cap that cannot be transferred back to the surface during the spark development 
time scale (see Appendix B in GM02) . The heat flow coefficient k influenced the required sur- 
face magnetic field B s in proportion to k ' 57 (for CR case). GM02 argued that for T s ~ 1 MK 
this effect could reduce the required B s by a factor of about 2. However, if T s ~ 3 MK as 
indicated by the case of PSR B0943+10, the reduction effect is negligible and we ignored it 
in this paper. 

8. Discussion and Conclusions 

The phenomenon of drifting subpulses has been widely regarded as a powerful tool 
for understanding the mechanism of coherent pulsar radio emission. RS75 first proposed 
that drifting subpulses are related to E x B drifting sparks discharging the high potential 
drop within the inner acceleration region above the polar cap. The subpulse-associated 
streams of secondary electron-positron plasma created by sparks were penetrated by much 
faster primary beam. This system was supposed to undergo a two-stream instability, which 
should lead to generation of the coherent radiation at radio wave lengths. However, careful 
calculations of the binding energy (critical for spark ignition) and the growth rate of the 
two-stream instability have shown that neither the sparking discharge nor the two-stream 
instability were able to work in a way proposed by RS75. Nonetheless, qualitatively their idea 
was still considered attractive, at least to the authors of this paper who have been continuing 
efforts to search for mechanisms that would actually make the RS75 model to work. In this 
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paper we applied the partially screened gap model proposed by GMG03 to PSR B0943+10, 
a famous drifter for which ZSP06 successfully attempted to measure the thermal X-ray from 
hot polar cap. We derived a simple relationship between the X-ray luminosity L x from the 
polar cap heated by sparks and the tertiary periodicity P3 of the spark-associated subpulse 
drift observed in radio band. This relationship reflects the fact that both the drifting (radio) 
and the heating (X-rays) rates are determined by the same value of the electric field in 
the partially screened gap. As a consequence of this coupling equations (13) and (14) are 
independent of details of the acceleration region. In PSRs B0943+10 and B1133+16, the 
only two pulsars for which both L x and P3 are measured, the predicted relationship between 
observational quantities holds very well. We note that P 3 is very difficult to measure and 
its value is known only for four pulsars: B0943+10, B1133+16, B0826-34 and B0834+06. 
The successful confrontation of the predicted X-ray luminosity with the observations in 
PSRs B0943+10 and B1133+16 encourages further tests of the model with future X-ray 
observations of other drifting pulsars. This is particularly relevant to PSR B0834+06, whose 
predicted X-ray luminosity is much higher than in PSR B0943+10, while the distance to 
both pulsars is almost identical. It is worth mentioning that due to relatively poor photon 
statistics it is still not absolutely clear whether the X-ray radiation associated with polar cap 
of PSR B0943+10 is thermal or magnetospheric (or both) in origin. However, the case of 
Geminga pulsar B0633+17 strongly supports thermal origin. Observations of PSR B0834+06 
with XMM-Newton should help to resolve this question ultimately. 

In both the steady SCLF model and the pure vacuum gap model, the potential increases 
with height quadratically at lower altitudes. However the growth rate is significantly different 
- the latter is faster by a factor of R/r p . It is well known that in the pure vacuum case (RS75), 
the potential grows so fast that a primary particle could quickly generate pairs with a high 
multiplicity, and that some backward returning electrons generate more pairs and soon a 
"pair avalanche" occurs and the potential is short out by a spark. In the PSG model we 
are advocating, the potential drops by a factor of r\. For essentially all the cases we are 
discussing, this rj value is much larger than the r p /R value required for the steady SCLF 
to operate, although it is much less than unity. The steady state condition is not satisfied, 
and the the gap is more analogous to a vacuum gap, i.e. the pair discharging happening 
intermittently. In fact, the partially screened potential drop is still above the threshold 
for the magnetic pair production, which in strong and curved surface magnetic field is a 
condition necessary and sufficient for the sparking breakdown. 

The original RS75 pure VG model predicts much too high a subpulse drift rate and 
an X-ray luminosity to explain the case of PSR B0943+10 and other similar cases. Other 
available acceleration models predict too low a luminosity and the explanation of drifting 
subpulse phenomenon is generally not clear at all (see ZSP05 for more detailed discussion). 
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In summary, the bolometric X-ray luminosity for the space charge limited flow (Arons & 
Scharlemann 1979) pure vacuum gap (RS75) and partially screened gap (GMG03) is (10~ 4 — 
10~ 5 )E (Harding & Muslimov 2002), (10" 1 - 10~ 2 )E (ZSP05), and ~ 10~ 3 E (this paper), 
respectively. The latter model also predicts right E x B plasma drift rate. Thus, combined 
radio and X-ray data are consistent only with the partially screened VG model, which requires 
very strong (generally non-dipolar) surface magnetic fields. Observations of the hot-spot 
thermal radiation almost always indicate bolometric polar cap radius much smaller than the 
canonical Goldreich- Julian value. Most probably such a significant reduction of the polar cap 
size is caused by the flux conservation of the non-dipolar surface magnetic fields connecting 
with the open dipolar magnetic field lines at distances much larger than the neutron star 
radius. 

Our analysis suggests the following pulsar picture: In the strong magnetic fields ap- 
proaching 10 14 G at the neutron star surface, the binding energy is high enough to prevent 
a full thermionic flow from the hot polar cap at the corotation limited level. A partially 
screened vacuum gap develops with the acceleration potential drop exceeding the threshold 
for the magnetic pair formation. The growth of this potential drop should be naturally lim- 
ited by a number of isolated electron-positron spark discharges. As a consequence, the polar 
cap surface is heated by back-flowing particles to temperatures T s ~ 10 6 K, just below the 
critical temperature T c at which the thermionic flow screens the gap completely. The typical 
radii of curvature of the field lines 1Z is of the order of polar cap radii r p ~ 10 3 — 10 4 cm. The 
only parameter that is thermostatically adjusted in a given pulsar is the shielding parameter 
7] = 10-3(^/^)^6 a5 P ~ 0.001T 6 L43 ftg a5 P, which determines the actual level of charge 
depletion with respect to the pure vacuum case (77=1) , and in consequence the polar cap 
heating rate as well as the spark drifting rate. It is worth to emphasize that 77 ~ 0.1 for 
longer period pulsars ?~ 1 s and rj ~ 0.01 for shorter period pulsars. Our calculations are 
consistent with PSR B0943+10 and few other drifting pulsars, for which the signatures of 
X-ray emission from the hot polar cap were detected. 

The sparks operating at the polar cap generate streams of secondary electron-positron 
plasma flowing through the magnetosphere. These streams are likely to generate beams of 
coherent radio emission that can be observed in the form of subpulses. Usov (1987) first 
pointed out that the non-stationarity associated with sparking discharges naturally leads to 
a two-stream instability as the result of mutual penetration between the slower and the faster 
plasma components. Asseo & Melikidze (1998) developed this idea further, calculated the 
growth rates, and demonstrated that the instability is very efficient in generating Langmuir 
plasma waves at distances of many stellar radii r ins ~ \0 A ~ 5 h CR = 10 7 ~ 8 cm, where h CR is 
the height of the acceleration region described by equation (21) (Melikidze, Gil & Pataraya 
2000). This is exactly where pulsar's radio emission is supposed to originate (e.g. Kijak, & Gil 
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1998). Conversion of these waves into coherent electromagnetic radiation escaping from the 
pulsar magnetosphere was considered and discussed by Melikidze, Gil & Pataraya (2000) 
and Gil, Lyubarski, & Melikidze (2004). These authors demonstrated that the nonlinear 
evolution of Langmuir oscillations developing in pulsar's magnetosphere leads to formation of 
charged, relativistic solitons, able to emit coherent curvature radiation at radio wavelengths. 
The component of this radiation that is polarized perpendicularly to the planes of dipolar 
magnetic field can escape from the magnetosphere (see Lai et al. 2001, for observational 
evidence of such polarization of pulsar radiation). The observed pulsar radiation in this 
picture is an indirect consequence of sparking discharges within the inner acceleration region 
just above the polar cap. In light of this paper we can therefore argue that the coherent 
pulsar radio emission is conditional on the presence of strong non-dipolar surface magnetic 
fields at the polar cap, with a strength about 10 13-14 G and radius of curvature of the order 
10 4 cm. 

In the very strong surface magnetic field assumed within the accelerator, processes such 
as photon splitting (Baring & Harding 2001) and bound pair creation (Usov & Melrose 1995) 
may become important. It has been suggested that such processes could potentially delay 
pair creation and thus increase the height and voltage of the accelerator. For the photon 
splitting case, the delay is significant only if photons with both polarization modes split - a 
hypothesis in strong magnetic fields (Baring & Harding 2001). It could be possible that only 
one mode split (Usov 2002). In such a case the gap height and potential of a PSG would not 
be significantly affected due to the high pair multiplicity in strong, curved magnetic fields. 
For bound pairs (Usov & Melrose 1995, 1996), in the very hot environment near the neutron 
star surface (with temperatures exceeding 1 MK), it is possible that bound pairs could not 
survive long from photoionization. Following Bhatia, Chopra, & Panchapakesan (1992) and 
Usov & Melrose (1996) one can roughly estimate the mean free path for bound pair photo- 
ionization. It turns out that for temperatures around (2-3) MK this mean free path is of 
the order of few meters, which is considerably less than the height of CR-driven PSG. In 
such a case, even if the bound pairs are initially produced, they would not significantly delay 
the pair formation. One could address this potential potential problem by referring to the 
detailed case study of PSR B0943+10 analyzed in this paper. This is the only pulsar in 
which we have full information concerning E x B drift rate and the polar cap heating rate. 
In the analysis we have used two methods. The first method is independent on details of 
accelerating region (such as height, potential drop, etc.) but is based only on the subpulse 
drift radio-data (eqs. [13,14]). The second one includes the detailed treatments of model 
parameters without considering the delaying effect by photon splitting and bound pairs 
(eqs. [25,26]). Both methods give consistent results, as illustrated in Figures (1-4). We 
therefore conclude that the delaying processes, if any, are not significant in the PSGs given 
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the physical conditions we invoke. Finally, we note that Usov & Melrose (1995) presented a 
model of the inner accelerator that shares some similar features with our model, although it 
is basically different. In their model the potential drop is self-regulated by partial screening 
close to the threshold for field ionization of bound pairs, so the surface temperature is 
maintained at the level necessary for this screening. In our model, the partial screening keeps 
the surface temperature slightly below the critical temperature at which the thermionic flow 
is space charge limited Arons & Scharlemann (1979). In the light of results presented in this 
paper we claim that bound pairs do not affect the formation of inner accelerators in pulsars. 

We acknowledge the support of the Polish State Committee for scientific research under 
Grant 2 P03D 029 26 (JG & GM) and NASA grants NNG05GB67G, NNG05GH91G and 
NNG05GH92G (BZ). 
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Fig. 1. — Models of partially screened inner acceleration regions (PSG) driven by curvature 
radiation seed photons above positively charged (ion case) polar cap for 102 pulsars from 
Table 1, sorted according to pulsar period. The horizontal axes correspond to the pulsar 
number Nl (bottom) or pulsar period (top). The vertical axes correspond to the surface 
magnetic field B s / B q (left-hand side), or surface temperature T s /10 6 K (right-hand side). 
The calculations were made for conditions corresponding to very strong (B s > 5 x 10 12 G) 
and curved (0.1 > TZq > 0.005) surface magnetic field. The vertical lines correspond to the 
case of PSR B0943+10 (Nl=41) and the hatched area encompassing three horizontal lines 
correspond to the range of surface magnetic field and temperature inferred for this pulsar 
from the XMM-Newton observation (ZSP05). These models allow to read off the physical 
conditions existing in the acceleration region above the polar cap of a particular pulsar in 
the form of parameters such as B s , TZq and rj. 
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Fig. 2. — Dependence of the surface temperature on the surface magnetic field in the case of 
NTPSG for CR (upper panel) and ICS (lower panel) cases. Solid lines represent the critical 
ion temperature 7], while the symbol-broken lines represent the actual surface temperature 
T s for different combinations of the parameters TZq and rj. As one can see only the CR case 
can work in pulsars with positively charged (ion case) polar cap. The solution for PSR 
B0943+10 marked by the dashed lines is TL % = 0.009, r) = 0.09 and T s = 3.1 MK, implying 
B s = 4.8B q ~ 2 x 10 14 G. 
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Fig. 3. — The thermal radiation efficiency from the hot polar cap for 102 pulsars from 
Table 1 for different parameters of the PSG marked in the legends. Pulsars are sorted 
according to the pulsar spin-down luminosity (N2 in Table 1). For clarity of presentation the 
surface field is fixed at B s / B q =4.8. The vertical dashed lines (N2=56) and three horizontal 
dashed lines (L x = 5 ± 2 x 10 28 erg/s) within the hatched area correspond to the special 
case of PSR B0943+10 (Tables 1 and 2). Three other pulsars for which values of P 3 /P is 
directly determined from single pulse radio data are included as diamond symbols: B0826-34 
(N2=10), B1133+16 (N2=51) and B0834+06 (N2=60). 
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Fig. 4. — Surface temperature of hot polar cap for 102 pulsars from Table 1 for different 
parameters of the PSG marked in the legends. Pulsars are sorted according to the pulsar 
period (Nl in Table 1). The vertical dashed lines (Nl=41) and three horizontal dashed lines 
(T s = 3 ± 1 x 10 6 K) within the hatched area correspond to the case of PSR B0943+10. 
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Table 1: Comparison of observed and predicted parameters 
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Table 2: List of 102 pulsars with drifting subpulses compiled from R86, GM02 and WES05. 
Nl and N2 refer to sorting according to pulsar period P and spin-down luminosity E, re- 
spectively. 
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